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ABSTRACT 

We describe a new spectroscopic technique for measuring radial metallicity gradients 
out to large galactocentric radii. We use the DEIMOS multi-object spectrograph on 
the Keck telescope and the galaxy spectrum extraction technique of Proctor et al. 
(2009). We also make use of the metallicity sensitive near-infrared (NIR) Calcium ii 
triplet (CaT) features together with single stellar population models to obtain metal- 
licities. Our technique is applied as a pilot study to a sample of three relatively nearby 
(< 30 Mpc) intermediate-mass to massive early-type galaxies. Results are compared 
with previous literature inner region values and generally show good agreement. We 
also include a comparison with profiles from dissipational disk-disk major merger 
simulations. Based on our new extended metallicity gradients combined with other 
observational evidence and theoretical predictions, we discuss possible formation sce- 
narios for the galaxies in our sample. The limitations of our new technique are also 
discussed. 
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' 1 INTRODUCTION 



The study of the stellar population parameters of early-type 
galaxies (ETGs) is an invaluable tool to understand the star 
formation and enrichment history of such galaxies. Unfor- 
tunately, the determination of these parameters is greatly 
hindered by the age-metallicity degeneracy (Worthey 1994) . 
This degeneracy affects photometric studies to a higher de- 
gree than spectroscopic studies. Because it is observationally 
expensive, typical spectroscopic studies of stellar popula- 
tions concentrate on the central regions of galaxies reaching 
out to < Ire, where is the effective radius. This implies 
that typically at least 50 per cent of a galaxy's stellar mass 
is not surveyed spectroscopically, leaving an incomplete pic- 
ture. 

Moreover, feedback processes such as galactic winds 
from supernovae (SNe) or stellar mass loss, and active galac- 
tic nuclei (AGN) play an important role in the star formation 
history of galaxies. This is because feedback processes eject 
the gas that is required for the birth of new stars thereby 
shutting down star formation. Because metals are produced 
in stars, the shut down of star formation and thus feedback 
processes have a direct influence on the enrichment history 
of a galaxy. The efflciency of those feedback processes varies 



with galactocentric radius. On the one hand, galactic winds 
are more efficient at expelling gas at large galactocentric 
radii where the local potential well is shallower (Matteucci 
1994; MartineUi, Matteucci & Colafrancesco 1998). On the 
other other hand, due to their location and high efficiency, 
AGN can quench star formation at small and possibly out 
to large radii (e.g., Croton et al. 2006; Wang & Kauffmann 
2008). Therefore, radial profiles of stellar population param- 
eters such as metallicity gradients can differ greatly from one 
galaxy to another depending on the relative importance of 
these feedback mechanisms and the formation mode (e.g., 
Spolaor et al. 2009). 

The processes involved in the formation and evolution 
of galaxies are still largely unknown despite tremendous ef- 
forts. The most popular scenarios have been the monolithic 
collapse and the hierarchical merging scenarios. In the classi- 
cal monolithic collapse scenario (Larson 1974; Larson 1975; 
Carlberg 1984; Arimoto & Yoshii 1987), galaxies form early 
via the dissipation of large gas clouds, which induces the 
formation of the vast majority of the stars. Enriched gas 
from rapidly evolving stars in the outskirts sinks towards 
the centre thereby causing steep radial metallicity gradients 
(i.e., decreasing metallicity with galactocentric radius). The 
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initial star formation period is then followed by quiescent 
evolution punctuated by few possible mergers and little sub- 
sequent star formation. 

In contrast, in the classical hierarchical scenario (Kauff- 
mann, White & Guiderdoni 1993; Kauffmann & Chariot 
1998; De Lucia et al. 2006; De Lucia & Blaizot 2007) galaxies 
form and evolve via the merging of sub-units. Star forma- 
tion is then triggered by merger events throughout a galaxy's 
history. Under this paradigm, the assembly of galaxies is a 
continuous process whereas the scale of the merger induced 
star formation episodes itself depends on the gas content of 
the progenitors involved. The high merger rates predicted by 
the hierarchical scenario should weaken or wash out metal- 
licity gradients (White 1980, but see Di Matteo et al. 2009). 
However, as recently shown by the simulations of Hopkins 
et al. (2008, hereafter H08), merger induced central star for- 
mation due to dissipation in gas-rich mergers can also create 
both shallow and steep radial metallicity gradients (see also 
Kobayashi 2004). 

Over the years these two scenarios have been refined 
and modified to accommodate observations and their pre- 
dictions accordingly updated. Recent literature has shown 
that modern versions of each of these scenarios can success- 
fully reproduce elliptical-like objects (Kampakoglou, Trotta 
& Silk 2008) . At large radii the shape of the radial metallic- 
ity gradient predicted by these scenarios differs. The effect 
of dissipation, such as in a gas rich merger or under the 
monolithic collapse paradigm, yields steep gradients in the 
central regions of galaxies. However, at large radii the pre- 
dictions for the slope of the metallicity gradient in gas-rich 
merger simulations tend to flatten (H08) while in the mono- 
lithic collapse simulations significant drops in metallicities 
are sometimes present (see Pipino, D'Ercole & Matteucci 

2008) . In the absence of dissipation, such as in dry merg- 
ers, the predicted gradients are shallower and more or less 
constant with radius. Therefore, the change in slope of the 
metallicity gradient contains a 'footprint' of the importance 
of dissipation in the star formation history of a galaxy. 

The techniques used to derive metallicity gradients of- 
ten suffer from the age- metallicity degeneracy. In their early 
study, Peletier & Valentijn (1989) used radial colour gra- 
dients as probes for metallicity gradients assuming a con- 
stant age. Colour-magnitude diagrams of resolved individual 
stars have also been used to more robustly determine the 
metallicity of halo star populations for very nearby galax- 
ies (e.g., Elson 1997; Harris & Harris 2002; Rejkuba et al. 
2005; Harris et al. 2007; Vlajic, Bland-Hawthorn & Freeman 

2009) . For unresolved stellar populations, even spectroscopic 
studies are not immune to the age-metallicity degeneracy as 
many authors have relied on only one or a few metal-lines 
measured from blue spectra (e.g., Spinrad et al. 1971; Cohen 
1979; Gorgas, Efstathiou & Aragon Salamanca 1990; Car- 
oUo, Danziger & Buson 1993; Davies, Sadler & Peletier 1993; 
Kobayashi & Arimoto 1999; Ogando et al. 2005; Forbes, 
Sanchez-Blazquez & Proctor 2005; Weijmans et al. 2009). 
In order to fully break the age-metallicity degeneracy, sev- 
eral studies (e.g., Sanchez-Blazquez et al. 2007; Brough et 
al. 2007; Spolaor et al. 2008b hereafter S08b) have used the 
technique of Proctor & Sansom (2002), which uses fitting 
with as many spectral indices as possible, thereby yielding 
more accurate metallicity values. Unfortunately, even with 
this technique both the finite slit length and the faintness 



of a galaxy at large galactocentric radii prevented previ- 
ous authors from efficiently measuring metallicity gradients 
to large radii especially in giant galaxies. Such studies have 
been recently compiled and expanded in Spolaor et al. (2009) 
to reveal an interesting inverse 'V-shape' relationship be- 
tween the mass of galaxies and the slope of their inner metal- 
licity gradient. At the low-mass end, the metallicity gradient 
goes from slightly positive slopes to more negative slopes as 
the mass increases. The relationship then reaches a turn- 
around at a mass corresponding to a velocity dispersion of 
roughly 140 km s~^ after which galaxies with higher velocity 
dispersion (or mass) tend to have less negative metallicity 
gradients. The brightest cluster galaxies do not seem to fit 
this trend however and can scatter from null down to large 
negative gradient values. 

In this pilot study, we extend the technique developed 
by Proctor et al. (2009, hereafter P09) to measure metallicity 
gradients out to large galactocentric radii for 3 giant ETGs 
(NGC 1407, NGC 2768 and NGC 4494). Following P09, we 
extract near-infrared (NIR) spectra of the galaxy halo light 
at large galactocentric radii using the DEIMOS multi-object 
spectrograph on Keck. The DEIMOS spectrograph is most 
efficient in the NIR where the Calcium ll triplet (CaT) spec- 
tral feature dominates (~ 8600 A). The CaT has been shown 
to correlate with metallicity (e.g., Armandroff & Zinn 1988; 
Diaz, Terlevich & Terlevich 1989; Cenarro et al. 2001) with 
little age sensitivity (e.g., Diaz, Terlevich & Terlevich 1989; 
Schiavon, Barbuy & Bruzual 2000; Cole et al. 2004; Carrera 
et al. 2007; Vazdekis et al. 2007, hereafter V03). It is possible 
that different Here we use the CaT together with the single 
stellar population (SSP) models of V03 to develop a new 
technique for deriving metallicity gradients using DEIMOS 
reaching as far out as ~ 2re. This technique can be applied 
to a large sample of galaxies and compared to the predictions 
of galaxy formation scenarios. 

Our paper is divided as follows: in Section [5] and [3] 
we give a brief description of our sample galaxies and an 
overview of our data, respectively. In Section |4] we explain 
the method used to extract metallicities and in Section[5]we 
give our results. Finally, Sections [6] and [7] contain a discus- 
sion of our results and a summary of this work, respectively. 



2 SAMPLE 

A summary of the relevant properties of our sample galaxies 
is given in Table [T] Digitized Sky Survey (DSS) images are 
shown in Fig. [1] Below we give a brief overview of each 
galaxy. 

2.1 NGC 1407 

NGC 1407 is the brightest group galaxy that dominates the 
Eridanus A group (Brough et al. 2006). It is a giant elliptical 
galaxy with a clear core-like luminosity profile (Spolaor et 
al. 2008a). It has been measured to have a uniform old age 
with a steep metallicity gradient and high a-element ratio 
within one effective radius (S08b). In their study of the stel- 
lar kinematics, Spolaor et al. (2008a) found tentative signs 
of a kinematically decoupled core (KDC) at the centre of 
NGC 1407. However, as they point out this KDC detection 
may be spurious and due to a possible misalignment of the 
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Galaxy 


Hubble 


P.A. 


Axis ratio 


Distance 


re 


Mb 


Mk 


Stellar mass 








Type 


(deg) 


(K band) 


(Mpc) 


(arcsec) 


(mag) 


(mag) 


(10" Mq) 


(km s"-"-) 


(km s-i) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC 1407 


EO 


60 


0.95 


26.8 


70 


-21.4 


-25.4 


2.86 


1782 


272.5 


NGC 2768 


SOi/2 


93 


0.46 


20.8 


64 


-20.9 


-24.6 


1.37 


1327 


181.8 


NGC 4494 


El-2 


173 


0.87 


15.8 


49 


-20.4 


-24.8 


1.64 


1335 


150.2 



Table 1. Galaxy properties. Hubble types (column 2) are as per NASA/IPAC Extragalactic Database (NED). Position angles and axis 
ratios (columns 3, 4) are from 2MASS. Distances (column 5) are based on surface brightness fluctuations (Tonry et al. 2001) and include 
the distance moduli correction of Jensen et al (2003). Effective radii (column 6) are taken from the Third Reference Catalogue of Bright 
Galaxies (RC3, de Vaucouleurs et al. 1991). B- and ii'-band absolute magnitudes (columns 7, 8) are calculated from RC3 and 2MASS 
apparent magnitudes, respectively, and using the distances quoted in column 5. Stellar masses (column 9) are calculated from the K-hand 
magnitude of column 8 assuming a M/Lk ratio corresponding to the V03 SSP of age 10 Gyrs and solar metallicity. Systemic velocities 
(column 10) are from P09. Central velocity dispersions (column 11) are as per Paturel et al. (2003). 




Figure 1. Images of our sample galaxies from the DSS with over- 
laid position of the slits whose spectra have signal-to-noise ratio 
> 8.5 (small circles). Slits identified as crosses are discussed in 
Section l4.1l Large circles correspond to 1, 2 and 3 effective radii. 
North and East are towards the top and left, respectively. 



slit with the semi-major axis of the galaxy. Nevertheless, 
if confirmed, the presence of a KDC in the centre of NGC 
1407 may be the signature of a merger history. Otherwise, 
it shows little sign of fine structure or recent disturbance. 



2.2 NGC 2768 

The literature is not unanimous with regards to the morpho- 
logical classification of NGC 2768, which ranges from ellipti- 
cal (de Vaucouleurs et al. 1991) to lenticular (Sandage, Tam- 
mann & van den Bergh 1981; Sandage & Bedke 1994). Evi- 
dence of rotation has been found by several authors (Fried & 
Ilhngworth 1994; Emsellem et al. 2004; P09) as well as the 
presence of a centrally concentrated kinematic twist (Mc- 
Dermid et al. 2006). NGC 2768 is located at the centre of 
a loose group of which it is the brightest galaxy (Giuricin 
et al. 2000). No previous metallicity gradient measurement 
is available, however the measured central metallicity varies 
from [Fe/H] « -0.2 to +0.3 (see HoweU 2005; Denicolo 



et al. 2005; Sil'Chenko 2006). The literature is also dis- 
crepant with respect to its central age with ages ranging 
between 8 and 15 Gyrs (e.g. Howell 2005; Denicolo et al. 
2005; Sil'Chenko 2006) possibly because of different spatial 
sampling. 



2.3 NGC 4494 

NGC 4494 is classified as an elliptical galaxy located in the 
Coma I cloud (Forbes et al. 1996; Larsen et al. 2001). Its 
luminosity profile displays a central cusp (Lauer et al. 2007). 
It has a KDC (Bender et al. 1988), an intermediate central 
age of 6.7 Gyrs and a central metallicity of [Fe/H] ~ +0.03 
(Denicolo et al. 2005). For these reasons, it is considered 
a good candidate for a possible gas-rich merger remnant 
(H08). 



3 DATA 

3.1 Acquisition 

Our main observing program was to obtain spectra for glob- 
ular cluster systems around ETGs. A total of 5 ETGs have 
been observed so far (see P09 for details). From this initial 
sample, we selected 3 galaxies that had the greatest number 
of spectra and highest quality (signal-to-noise ratio). 

Spectra were obtained using the DEIMOS spectrograph 
on the Keck telescope during the nights of 2006 November 
19-21, 2007 November 12-14 and 2008 April 8. The seeing 
was good (typically ~ 0.7 arcsec). The 1200 1 mm~^ grating 
was used with a central wavelength of either 7500 A (2006 
Nov.) or 7800 A (2007 Nov., 2008 Apr.). In every case, the 
setup allowed for the coverage of the CaT region (~ 8400 — 
8900 A) with a resolution of AA ~ 1.5A for the 1" slit width. 
A total of 6, 2 and 3 multi-slit masks were observed for NGC 
1407, NGC 2768 and NGC 4494, respectively. The typical 
total exposure time on each mask for NGC 1407 and NGC 
2768 was 2 hours (4 x 30 minutes exposures) and 1.5 hours 
(3 X 30 minutes exposures) for NGC 4494. One mask for 
NGC 2768 was observed for 1 hour on two separate nights 
yielding two independent measurements. Fig. [T] shows our 
selected galaxies together with the positions of the selected 
slits (see Section [3.21 for the selection process). 
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Galaxy 


Central passband 


Continuum passbands 




(A) 


(A) 


NGC 1407 


8605.0 - 8695.5 


8526.0 - 8536.0 






8813.0 - 8822.0 


NGC 2768 


8605.0 - 8695.5 


8478.0 - 8489.0 






8813.0 - 8822.0 


NGC 4494 


8605.0 - 8695.5 


8478.0 - 8489.0 






8813.0 - 8822.0 



1000 



Table 2. Sky scaling index definitions. 



3.2 Reduction 

The DEIMOS data were reduced using the IDL spec2d data 
reduction pipehne written for the DEEP2 galaxy survey. The 
pipehne performs both the flat fielding using internal quartz 
flats and the wavelength calibration using the ArKrNeXe 
arc lamps. Residual fringing is negligible due to instrumental 
design (e.g., Faber et al. 2008; Wirth et al. 2004). In addition 
to the sky subtracted globular cluster spectra, the pipeline 
produces several outputs, among which are the background 
(or 'sky') spectra. 

We used the technique of P09 to extract galaxy halo 
light spectra out of the background spectra. The background 
spectra are essentially the sum of both sky and galaxy light. 
P09 use an appropriately scaled 'true' sky spectrum that is 
then subtracted from the background spectrum to extract 
the sky subtracted galaxy halo light spectrum. 

In this work, we compute the scaling factor on the raw 
background spectra as the excess flux in the region 8605.0- 
8695. SA above a linear continuum determined from two 
carefully selected sidebands that avoid both strong skylines 
and galaxy spectral features for each galaxy (see Table 
Fig.|2]shows an example background spectrum together with 
the definition of the sky scaling factor. As can be seen in 
Fig. H the recession velocity of NGC 2768 (and NGC 4494) 
causes the Ca3 line to be partly redshifted into the central 
passband of the sky index. This could introduce systematic 
errors when applying the sky subtraction. To test this, we 
use a slightly narrower central passband avoiding the Ca3 
feature and find no noticeable changes to our results. Thus, 
because small variations in instrument resolution and wave- 
length solution across the mask could cause the minima be- 
tween the individual skyline peaks to vary and introduce 
errors if the edge of the sideband is near a strong skyline we 
choose to use the same central band definition for all three 
galaxies in our sample. 

The 'true' sky is derived from the normalized sum of 
several background spectra at large galactocentric radii (6- 
7re). Even at 6-7re, there is still some galaxy background 
light in our sky estimate. However, using a de Vaucouleurs' 
(1953) luminosity profile, we estimate that the galaxy light 
in these outer sky spectra is at most 6 percent of the galaxy 
light in our science spectra. We use Monte Carlo meth- 
ods to evaluate the accuracy of the sky subtraction at 9 
signal-to-noise intervals ranging from 10 to 50 and find 
that the error in the final continuum level introduced by 
over/undersubtracting the sky is negligible compared to the 
noise (0.7 percent of the level of the noise for our lowest 
signal-to-noise ratio) if a linear continuum is assumed as 
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8500 



8600 8700 
wavelength (A) 



8800 



Figure 2. Example of a background spectrum (i.e., galaxy light 
+ sky) for NGC 2768. Labelled gray shaded regions highlight the 
continuum and central passbands of the sky scaling index defini- 
tion. The dashed line shows the continuum estimate derived for 
this particular spectrum (the sky subtracted spectrum is shown 
in the top panel of Fig. [3]l. Positions of galaxy individual CaT 
features are labeled. 



is typical in this spectral region. The method described in 
P09 is similar to that of Norris et al. (2008) and Proctor et 
al. (2008), which employed background spectra from Gem- 
ini/GMOS and Keck/LRIS, respectively. For further infor- 
mation on the sky subtraction method used herein see P09. 

Next, the galaxy halo light spectra were fitted using 
the pPXF code of Capellari & Emsellem (2004) to extract 
halo kinematics out to large galactocentric radii (< Sr^). 
The kinematics are presented in P09. In Fig. O we show ex- 
amples of the extracted galaxy halo light spectra and high- 
light regions that are still contaminated by skyline residuals. 
We find that approximately 10 per cent of the amplitude of 
these skyline residuals can be attributed to small variations 
of the wavelength solution across the mask. The remain- 
ing fraction is likely caused by the inherent complications 
associated with non-local sky subtraction due to variations 
of the sky spectrum over time and across the large field of 
view. As noted in P09, the strong skyline residuals are not 
significantly larger than the Poisson noise associated with 
them. Because of these skyline residuals, our method yields 
better results for galaxies with systemic recession velocities 
that do not cause the CaT features to be shifted into sky- 
line dominated regions (i.e., 500 km s"^ < Vsys < 1400 km 
s"^ or 1700 km s"^ < Vsys < 2500 km s"^). Due to the in- 
herent difficulties associated with flux calibrating multi-slit 
data our spectra are not flux calibrated. 

Finally, we select our highest signal-to-noise spectra by 
removing any spectra with an average number of counts per 
angstrom < 73. This roughly corresponds to a signal-to- 
noise cut of 8.5. 



4 ANALYSIS 

4.1 Index Measurements 

There are several CaT index deflnitions that apply to inte- 
grated light spectra. The choice of index deflnition is some- 
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8500 8600 8700 



wavelength (A) 

Figure 3. Examples of a high (top panel), typical (middle panel) 
and low (bottom panel) signal-to-noise galaxy halo light spectra 
for NGC 2768. The raw spectra are shown as the black solid 
line. Gray highlighted regions show spectral ranges contaminated 
by skylines. Selected spectra have average number of counts per 
angstrom > 73. 



what dependent on the purpose one desires to fulfill. Indeed, 
ArmandrofI & Zinn (1988, hereafter AZ88) used narrow cen- 
tral passbands for the CaT to determine an empirical con- 
version between their CaT index and metallicities using the 
integrated light spectra of Galactic globular clusters. A year 
later, Diaz, Terlevich & Terlevich (1989, hereafter DTT89) 
defined a much broader index, which they applied to the in- 
tegrated light spectra of galaxies for which line indices are 
broadened due to their large velocity dispersion. The major 
drawback of this definition is that it uses the same con- 
tinuum passbands for all three CaT lines and is therefore 
very sensitive to variations in the shape of the continuum. 
In the work of Cenarro et al. (2001, hereafter COl), this is 
circumvented by the possibility of using an arbitrary num- 
ber of continuum passbands. Also, in order to reduce the 
effects of skylines and other non-Poisson noise during the 
CaT index measurement, the method of COl weighs each 
pixel according to its variance. Moreover, the three features 
that constitute the CaT (i.e., Cal, Ca2 and Ca3) have differ- 
ent relative depths with the bluest (Cal) being the weakest. 
For this reason, some authors have decided to give varying 
weights to the different CaT features or to remove the Cal 
feature altogether from their CaT index definition (e.g. Ar- 
mandroff & da Costa 1991; Rutledge et al. 1997; Koch et al. 
2006; Koch et al. 2008). 

Because we are using velocity dispersion broadened 
galaxy spectra, we cannot use the narrow AZ88 central 
passbands definition or their empirical conversion. We thus 



choose to adopt the DTT89 index central passbands def- 
inition with the continuum determination and index mea- 
surement technique described in COl. The central passband 
definitions as well as the continuum passbands of our CaT 
index are given in Table |3] and Fig. U We also employ a 
weighted sum such that: 

CaT = 0.4 X Cal + Ca2 + Ca3. (1) 

This is done in order to minimize the impact of the more 
uncertain Cal feature. The continuum passbands are chosen 
to uniformly cover the spectral range around the CaT fea- 
tures and in such a way as to avoid large spectral features 
and regions dominated by skyline residuals. The weights on 
the individual CaT features in Eq. [T]are chosen to minimize 
the estimated errors on the CaT index values. Error esti- 
mates are computed using the background spectra (i.e., be- 
fore sky subtraction) instead of fully propagated variance ar- 
rays, which are unavailable as our spectra are not processed 
within the data reduction pipeline. The background spectra 
provide a good estimate of the true variance since the bulk 
of the variance is caused by skylines, which in turn yields a 
robust error estimate. The technique described in Cardiel et 
al. (1998) and COl (appendix A2) for generic indices is used. 
Finally, we apply a velocity dispersion correction in order to 
obtain CaT index values that are comparable to those mea- 
sured at the models' dispersion. This correction is shown 
in Fig. [3 For both NGC 2768 and NGC 4494, the veloc- 
ity dispersions were taken from P09. However, the velocity 
dispersion profile of NGC 1407 measured by P09 reveals a 
sharp 'spike' around ~ 0.9re. The high velocity dispersion 
values measured from the corresponding slits (highlighted in 
white in Fig. [TJ produce high CaT index values. These in 
turn yield highly deviant (unphysical) metallicities of up to 
[Fe/H] ~ 2.0 around r ~ 0.9re. The cause of this 'spike' in 
the velocity dispersion profile is unknown. Because the re- 
quired velocity dispersion correction is clearly too large, we 
adopt velocity dispersion values extrapolated through the 
'spike' by assuming a smoothly declining velocity dispersion 
profile as typical of elliptical galaxies. 

The errors on the CaT index values are propagated ac- 
cordingly and include errors in the respective velocity dis- 
persion measurements. Because of the non-linear nature of 
the velocity dispersion correction, large velocity dispersions 
yield larger CaT index errors. Finally, we estimate that sky 
subtraction random errors contribute at most 0.001 A, which 
is minimal compared to our measured typical random errors. 
Larger systematics may arise from the skyline residuals as 
described in Section [3]2] These are alleviated and partly ac- 
counted for via weighting according to the COl index mea- 
surement method. In what follows we will refer to the CaT 
index measured using the above method as CaT. 

In order to facilitate comparison with previous studies, 
we use the V03 SSPs with 13 Gyrs to derive a conversion be- 
tween our CaT and the similar CaT index defined by DTT89 
(CaTnTT)- We obtain the following: 

CaToTT = 18.0 - 7.0 x CaT + 0.86 x CaT^ (2) 
with a standard deviation of cr = 0.25 A. 
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Figure 4. Representation of the CaT definition used in this worlc. 
The sohd line spectrum is a V03 SSP model with [Fe/H] = -0.68 
and age of 13 Gyrs broadened to 150 km s~^ dispersion to mimic 
the effect of a typical galaxy's velocity dispersion. Labelled gray 
shaded regions represent the CaT and continuum passbands. 
The dashed line shows the continuum level estimated using the 
method of COl. 



1.4 



1 .2 




Figure 5. Applied velocity dispersion index correction. The cor- 
rection was obtained by comparing our index measured directly 
on the V03 SSP spectra of ages 13 Gyrs (CaT) and on the same 
spectra convolved with Gaussians of corresponding widths for a 
given velocity dispersion {CaT„). The shaded area shows the ex- 
pected scatter in the correction caused by the different model 
metallicities. The solid line shows the best fit polynomial of third 
order through the median metallicity (i.e., [Fe/H] ~ —0.28) 
model points. 



Feature 


Weight 


Central passband 

(A) 


Shared continuum 
passbands (A) 


Cal 


0.4 


8483.0-8513.0 


8474.0-8483.0 


Ca2 


1.0 


8527.0-8557.0 


8514.0-8526.0 


Ca3 


1.0 


8647.0-8677.0 


8563.0-8577.0 








8619.0-8642.0 








8680.0-8705.0 



Table 3. CaT index definition. 



4.2 Conversion into metallicities 

We use the single stellar population (SSP) models of 
Vazdekis et al. (2003, hereafter V03) with the Kroupa (2001) 
initial mass function in order to convert our CaT measure- 
ments into metallicities. There are two main reasons that 
motivated this choice. First of all, the V03 models provide 
spectral energy distributions (SEDs) at a resolution compa- 
rable to that of our data. Moreover, the V03 models show 
good agreement with Galactic globular cluster data (see fig. 
14 of V03). We assume a constant age of 13 Gyrs to con- 
vert CaT into metallicity. The predictions of the V03 models 
cover metallicities ranging from —1.68 < [Fe/H] < -1-0.20 
(corresponding to 4.1 A < CaT < 6.2 A). Outside this 
range, the behavior of the CaT with metallicity is not con- 
strained. This can be a serious issue particularly at high 
values of CaT and [Fe/H] where the conversion is steepest. 
For this reason and to avoid the introduction of uncertainties 
related to model extrapolation, we apply a hard boundary 
such that metallicities derived from data points that have 
measured CaT > 6. 2 A are discarded and assigned the max- 
imum metallicity (i.e., [Fe/H] — 0.2 dex). These points are 
clearly identified in the plots that follow and are not used in 
subsequent analyses. Errors are propagated from the CaT 
errors. Once again, the non-linearity of the derived conver- 
sion from CaT into metallicity yields larger error bars at 
high metallicities. We also make use of rolling averages with 
radius as they are more robust against random fluctuations. 

There could be other sources of uncertainty. First, there 
could be hot blue stars such as young, blue horizontal branch 
(BHB) or blue straggler (BS) stars contaminating the CaT 
with their predominant Paschen line features (see COl). Un- 
fortunately, varying BHB morphologies and blue stragglers 
are not modeled by V03 and one has to worry about possi- 
ble contamination by hot blue stars when converting CaT 
into metallicities as 3 of the features in the Paschen series of 
hydrogen overlap with the 3 CaT lines. Ifowever, we visually 
inspected our spectra and find that there is no indication for 
the presence of a Paschen line at 875lA where it should be 
most easily seen. 

Another source of uncertainty in our conversion into 
metallicities has to do with the effect of age on the CaT 
feature. Fortunately, the CaT is only minimally influenced 
by age effects for ages > 2.5 Gyrs (e.g., DTT89; Schiavon 
et al. 2000; V03; Cole et al. 2004; Carrera et al. 2007). This 
is in apparent contradiction with the predictions of the V03 
models shown in Fig. [S] On the other hand, it is likely that 
different CaT index deflnitions could yield different age de- 
pendencies as is often the case for other spectral features. 
Thus, it is conceivable that the apparent age trends seen in 
Fig. [6] may be influenced by our choice of index definition. 
Nevertheless, if we ignore the points corresponding to the 5 
Gyrs models with metallicities [Fe/H] = —1.68 and -1.28 
that V03 consider unreliable due to a lack of corresponding 
stars in the stellar library, the maximum error on the in- 
ferred metallicity induced by wrongly using a 13 Gyrs old 
SSP to estimate metallicity in a 5 Gyrs SSP is insignificant, 
as expected from observational studies. Therefore, because 
the inner parts of our galaxies are measured to be much older 
than 2.5 Gyrs, age effects should not strongly influence our 
inferred metallicities. 
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Figure 6. Single stellar population model predictions for 13, 9 
and 5 Gyrs (solid, dashed and dotted lines, respectively) from 
V03. The spectral resolution is AA ~ 1.5A. The 5 Gyr model 
data are uncertain for metallicities [Fe/_ff] = —1.28 and -1.68. 



5 RESULTS 

The spatial distribution of slits on the DEIMOS mask was 
optimized for the study of GCs. For this reason, the po- 
sition from where our halo light spectra are extracted are 
distributed in a random manner around a given galaxy. In 
order to present our results in a way that is comparable to 
previous measurements (i.e., along the major or minor axes), 
we transform our galactocentric radii values into their spher- 
ical equivalent if the slit lies along the semi-major axis. To 
do this, we first define the effective change in right ascension 
and declination as 



Aq = {a slit — Ci galaxy) C,Ofii galaxy, and 
A5 = 5 slit — Sgalaxy, 



(3) 



respectively; where we have implicitly used the small-angle 
approximation. The subscripts 'slit' and 'galaxy' correspond 
to the position of the slit and the photometric centre of the 
galaxy in equatorial coordinates, respectively. Next we use 
the photometric axis ratio ib/a), position angle </>, and ef- 
fective radius re from the literature (see Table [TJ to convert 
the 'true' distance to the centre of the galaxy into a pseudo 
major-axis distance. Or mathematically, we apply the fol- 
lowing formula: 



(r/re) : 



(A5')2 



(4) 



where Aa' = (A5)cos<^ + (AQ)sin<j!>, and A5' = 
— (Aq) cos (j> + (A5) sin (j). Assuming that there is no signifi- 
cant change in metallicity along isophotes, our results should 
be comparable to the literature's values for slits aligned with 
the semi- major axis. If this assumption is partly incorrect, 
it will introduce scatter in the metallicity at a given radius. 
In what follows, all galactocentric radii have been obtained 
with this method. 

We plot the rolling average using 8 to 10 points depend- 
ing on the size of the dataset in Fig. [T] |8]and|9] The rolling 
average gives a clearer indication of the actual shape of the 
metallicity gradient in our galaxies by eliminating the con- 
fusion caused by the scatter due to random errors on our 



data points. We find that the slope of the CaT and metal- 
licity gradients changes with galactocentric radius and thus 
cannot be appropriately described as linear in log-log space. 
This makes their objective quantification difficult. 

The quantification of the metallicity gradient in the in- 
ner regions of galaxies using the slope of a fitted line in 
log-log space has been useful in previous works to quantify 
the steepness of the gradient and constrain theoretical mod- 
els (e.g. Spolaor et al. 2009; Kobayashi & Arimoto 1999). 
However, we do not fit linear relations to our metallicity 
gradients as this does not quantify them properly. Instead, 
we compare the entire profile with theoretical expectations 
(see Section O. 



5.1 NGC 1407 

This galaxy is the only galaxy in our sample for which we 
have previous literature gradient values available to compare 
our results to. Unfortunately, it is not our best dataset in 
that the radial coverage is limited. Nevertheless, as seen in 
Fig. [71 the general agreement of our CaT measured metallic- 
ities with the metallicities measured from long-slit observa- 
tions by S08b is reasonable. Although the individual values 
agree with the SOSb data within the errors at similar galac- 
tocentric radii, the rolling average shows a slight offset in 
metallicity (~ 0.3 dex). The SOSb data were obtained from 
optical spectra and calculated using the well tested fit- 
ting method of Lick indices of Proctor & Sansom (2002). 
Considering that the present study may be sampling differ- 
ent stellar populations since 1) we sample a broader range of 
position angles, and 2) our wavelength range is much redder 
than SOSb, this broad agreement is an indication that the 
metallicities measured from the CaT are reliable. Moreover, 
while we estimate that errors due to the lack of flux calibra- 
tion of our spectra are negligible compared to the quoted 
errors, we cannot completely rule out the possibility that 
systematic offsets could be present. 

Both the CaT and derived metallicities reveal a steep 
gradient falling to [Fe/H\ ~ -l.S at log(r/re) ~ 0.15. Al- 
though the radial range covered here is limited, the CaT 
and metallicity gradients show no clear sign of leveling off. 



5.2 NGC 2768 

As can be seen in Fig.[Sl the data for NGC 2768 cover a more 
extensive radial range than our NGC 1407 data. They also 
contain more individual measurements than for any other 
galaxy in our sample and exhibit smaller scatter. The central 
metallicity value measured by Howell (2005) is consistent 
with our innermost data points. A visual inspection of Fig. 
[S] reveals a steep average metallicity gradient in the inner 
parts that steepens slightly before becoming shallower at 
log(r/re) ~ 0.1 around a metallicity of \Fe/ H] = —1.0. 

Although not inconsistent (i.e., still within the quoted 
errors) , a few outer slits seem to exhibit CaT values beyond 
the upper limit permitted by V03 (see Fig. [8]). Because of the 
uncertainties related to model extrapolation and especially 
given the steepness of the predicted relationship between 
CaT and [Fe/H] at high CaT values, we consider these 
metallicity data unreliable. 
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Figure 7. CaT (top panel) and metallicity (lower panel) gradients for NGC 1407. The stars are from the long slit data presented in 
S08b and the filled circles represent our data. Upper arrow error bars are used when the CaT is beyond the upper limit (thin dashed 
lines) of the V03 SSP models. Thick solid lines are rolling averages. Profiles from the dissipational merger models of H08 are shown as 
thick gray lines (labels as per Table l4t . 



5.3 NGC 4494 



Because shorter exposure times were used for NGC 4494, 
the number and spatial coverage of our measurements are 
much smaller than for the previous two galaxies. Despite the 
large scatter in our measured metallicity values, the agree- 
ment with the central estimate from Denicolo et al. (2005) 
is good. Interestingly, the data show no visible gradient in 
the raw CaT measurements (see Fig.|9} and the CaT values 
seem to scatter about the upper limit of the V03 models. 
Moreover, the metallicity "gradient" seen in Fig. [9] is ficti- 
tious and caused by our exclusion of the randomly scattered 
points above the CaT upper limit. Indeed, the majority of 
the measured metallicity points in NGC 4494 are very high 
(i.e., consistent with [Fe/H] > 0.20 dex) within the galacto- 
centric radius probed and an accurate metallicity gradient 
could not be measured. For this reason, we cannot comment 
on the presence or absence of a metallicity gradient in this 
galaxy based on the current dataset. 



It is clear that longer exposure times and more data 
are required to properly constrain the metallicity gradient 
in this galaxy at large galactocentric radii. Moreover, as for 
NGC 2768, the large CaT values found indicate that the 
method presented herein may not be accurate at metallic- 
ities near solar and particularly beyond [Fe/H] = -f0.20. 
Luckily, the negative radial metallicity gradients present in 
most galaxies mean that usually metallicities at large radii 
are sub-solar. 




Figure 8. CaT (top panel) and metallicity (lower panel) gradi- 
ents for NGC 2768. The star shows the central (r < re/8) metal- 
licity as measured by Howell (2005) and the circles represent our 
data. Hollow circles and upper arrow error bars are used when 
the CaT is beyond the upper limit (thin dashed lines) of the V03 
SSP models. Thick solid lines are rolling averages. Profiles from 
the dissipational merger models of H08 are shown as thick gray 
lines (labels as per Table |4]l. 



6 DISCUSSION 

In this pilot study, we develop a new technique to measure 
metallicity gradients at large galactocentric radii by expand- 
ing on the technique presented in P09. Using the CaT and 
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Figure 9. CaT (top panel) and metallicity (lower panel) gradi- 
ents for NGC 4494. The star shows the central (r < re/8) metal- 
licity as measured by Denicolo et al. (2005) and the circles repre- 
sent our data. Hollow circles and upper arrow error bars are used 
when the CaT is beyond the upper limit (thin dashed lines) of 
the V03 SSP models. Thick solid lines are rolling averages. The 
metallicity gradient shown in the lower panel is fictitious (see 
text). Profiles from the dissipational merger models of H08 are 
shown as thick gray lines (labels as per Table |4}. 
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Table 4. Relevant properties of H08 dissipative major merger 
simulations (column 1) selected based on their matching of the 
luminosity profiles of our sample galaxies (column 2) . The quoted 
stellar mass (column 3) is that of the remnant. Central velocity 
dispersions (column 4) are median line-of-sight values measured 
within 1 Tf, (as per column 5). fga3 (column 6) is the mass fraction 
of the merging disks in the form of cold gas just before the final 
merger. 



the SSP models of V03, we obtain spectroscopic metallic- 
ity measurements out to 1.4re, 1.9re and 1.5re for NGC 
1407, NGC 2768 and NGC 4494, respectively. We find that 
the metallicity gradients are not well described as a straight 
hne in log-log space as their slope varies with galactocentric 
radius. The slope variations themselves contain information 
about the processes involved in the galaxy's formation. 

In order to better understand our observed metallicity 
gradients and the processes involved in their formation, we 
compare our results to the predictions from the models of 
H08. In HOS, suites of hydrodynamic simulations of dissipa- 



tional disk-disk galaxy merger remnants are compared to the 
observed properties (sizes and surface brightness profiles) of 
individual elliptical galaxies. The simulations determine the 
enrichment (i.e. metallicity) self-consistently from star for- 
mation and include a prescription for AGN and supernova 
feedback. For the individual galaxies in our sample, the mod- 
els producing a remnant whose simulated surface brightness 
profile yield the best matches to a given galaxy's observed 
surface brightness profile are selected (see Table |3| with no 
prior on the metallicity or mass. Thus the predicted HOS 
metalhcity gradients are not scaled to match the observed 
data. In Fig. [T] [S] and [51 we show the line-of-sight averaged 
and B-band luminosity-weighted stellar metallicity profiles 
predicted 3 Gyrs after the merger of these selected models 
for NGC 1407, NGC 2768 and NGC 4494, respectively The 
simulated metallicity profiles are not strongly sensitive to 
the chosen line-of-sight or age of the merger. 

Although the model metallicities are determined self- 
consistently from stars formed in the merger simulations, 
progenitor disks at the start of the simulation must have 
some 'initial' metallicity. The chosen initial metallicities of 
the progenitors will not change the central predicted metal- 
licity of the remnant significantly as it is dominated by the 
merger induced dissipational star formation. However, the 
initial metallicity of the progenitors will come to dominate 
the metallicity of the remnant at sufficiently large radii. This 
sets an effective metallicity 'floor' at large radii. 

First we compare our results with the simulated sys- 
tems from HOS where the progenitor disks are initialized 
to lie on the observed redshift 2 = mass-metallicity rela- 
tion (e.g., Tremonti et al. 2004) with a uniform metallicity 
at all radii. We find that this predicts too high a metal- 
licity at large radii in the remnants to be compared with 
the current dataset. In the simulations shown in Fig. [T] [S] 
and [9l we lower the metallicity of the progenitors by a fac- 
tor of 3 (i.e., ~ —0.48 dex). Again the metallicities at 
< Ire are unchanged, but those at large radii are reduced. 
Under the dissipational major-merger formation paradigm, 
this lowered 'floor' suggests that 1) the remnant's metallic- 
ity increasingly reflects that of its progenitors' stars as the 
radius probed increases; and 2) since the galaxies in ques- 
tion have relatively early formation times {z > 1 to 2), the 
appropriate 'initial' mass-metallicity relation is not that ob- 
served today (i.e. z — 0) but that at these redshifts, which 
yields lower metallicities (e.g., Lara-Lopez et al. 2009). 

From the H08 simulations, we flnd that the central lu- 
minosity profile of NGC 1407 is well matched by mergers 
of intermediate gas-richness with gas fractions (fgas) at the 
time of merger ~ 10%, typical of L* galaxies. The theo- 
retical expectation is a relatively smooth metallicity profile, 
similar to that observed for log(r/re) < —0.1. Once again, 
the different profiles shown correspond to the three simula- 
tions that provide the best match to the observed surface 
brightness profile with no prior on the metallicity or mass 
of NGC 1407. Thus, the amount of dissipation needed to 
match the luminosity profile appears to also provide a good 
match to the metallicity gradients for log(r/re) < —0.1. Be- 
yond that radius however there is an observed steepening of 
the metallicity gradient, which is not predicted by the mod- 
els. This could be an indication that either 1) the initial 
metallicity of the progenitors is still too high (i.e., the 'fioor' 
should be set lower initially); or 2) that a dissipative ma- 
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jor merger remnant is not a good description of the forma- 
tion and evolution of NGC 1407. Under both the monolithic 
collapse scenario and dissipative merger scenarios, an early 
formation and assembly is required. Indeed, as concluded 
by S08b, NGC 1407's steep metallicity gradient, uniform 
old age and smooth photometric profile are consistent with 
a formation in which early dissipation played a major role. 
Moreover, it is not clear that a disk-disk major merger with 
plausible progenitors at moderate redshift could reproduce 
such a gradient unless the progenitors themselves possessed 
strong metallicity gradients (see Di Matteo et al. 2009). Al- 
ternatively, and still consistent with a hierarchical merging 
scenario, it is possible that the outskirts of NGC 1407 were 
built up primarily from smaller shredded systems (i.e., mul- 
tiple minor-mergers, see Naab, Johansson & Ostriker 2009), 
which typically have lower metallicities. 

The matching of the observed and simulated photomet- 
ric profiles of H08 is more ambiguous for NGC 2768 because 
it contains a lot of dust and exhibits several fine structures 
such as filaments and a ring of ionized gas (Martel et al. 
2004; Lauer et al. 2005). This causes ambiguity with re- 
spect to selecting the best fitting models. Indeed, H08 find 
reasonable matches with both relatively gas-poor mergers 
if gas ~ 3 — 5%) to relatively gas-rich (fgas 20 — 30%) ones. 
Unfortunately, the metallicity profiles are not sufficiently ac- 
curate to break this degeneracy. Nevertheless, the rolling av- 
erage suggests that the simulations of H08 are a reasonable 
match to the measured metallicity profile and particularly 
for the relatively gas-rich models (i.e., fgas ~ 20 — 30%). 

Unfortunately, the metallicity gradient for NGC 4494 
is not well constrained by our data. The shorter exposure 
time has yielded fewer data points to constrain the gra- 
dient. Moreover, at large CaT (and velocity dispersions), 
the inferred metallicities are more uncertain and may yield 
unreliable results especially in modest signal-to-noise data. 
For this reason, we could not detect and measure a reliable 
metallicity gradient in NGC 4494. On the other hand, our 
data suggest that the metallicity of NGC 4494 is very high 
(i.e., [Fe/H] > 0.20 dex) out to log{r/r^) ^ 0.2. 



7 SUMMARY AND FUTURE WORK 

We describe a new technique for obtaining radial metallicity 
gradients of galaxies out to large galactocentric radii using 
the DEIMOS multi-object spectrograph on the Keck tele- 
scope. We use the NIR CaT spectral feature and convert 
our CaT index into metallicity with the use of the V03 SSP 
models. This new technique is then applied to three inter- 
mediate mass to massive ETGs as a pilot study. Our results 
agree well with previous literature inner values. We find that 
at large galactocentric radii our measured metallicity gradi- 
ents are not well described with straight lines (in log-log 
space) and show significant variations with galactocentric 
radius. A comparison to theoretical models is used in order 
to interpret our metallicity gradients. 

We conclude that NGC 1407 is likely to have formed 
the bulk of its stars via dissipational processes. This is con- 
sistent with the monolithic collapse scenario. However, we 
also hypothesize that the low metallicity stellar populations 
probed in the outskirts of NGC 1407 could have been hier- 
archically assembled from smaller low-metallicity systems. 



For NGC 2768, the measured metallicity profile is well re- 
produced by the dissipative major merger models of H08. 
We were unable to measure the metallicity gradient in NGC 
4494 due to both the lower quality of this dataset and the 
fact that the measured CaT values scatter around the up- 
per limit permitted by the models. Nevertheless, we are able 
to conclude that the metallicity of NGC 4494 as measured 
via the CaT is [Fe/H] > 0.20 dex within the galactocen- 
tric radii probed. The limitations of our new technique as a 
function of signal-to-noise, metallicity and velocity disper- 
sion are also discussed. 

The main science goal of this dataset was to obtain 
spectra of GCs. In the future we plan to include dedicated 
galaxy halo background slits in order to increase the ra- 
dial and azimuthal coverage. This will enable 2-dimensional 
mapping of the metallicity distribution similar to the works 
of Kuntschner et al. (2006) and Rawle et al. (2008) reach- 
ing out to larger radii due to the larger field of view of the 
DEIMOS spectrograph. The assumption that the distribu- 
tion of the metallicity follows the isophotes, which was im- 
posed here, could then be tested and its effect on our con- 
clusions assessed. 
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APPENDIX A: DATA TABLES 

This section is a compilation of the results plotted in Fig. 
[3 ID and m We quote the maximum value of [Fe/H] = +0.2 
when the measured CaT is larger than the range covered 
by the V03 models. However, we wish to caution the reader 
that while the quoted CaT values are reliable within the 
quoted errors, the accuracy of the metallicity estimates is 
unknown for these specific data points only. 
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Table Al. Individual values for NGC 1407. Columns 1 and 2 give 
the position of the individual slits in right ascension and decli- 
nation (12000), respectively. The major axis equivalent radius is 
given in column 3. Measured values of the CaT index and [Fe/H] 
are shown in columns 4 and 5. 
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Table A2. Individual values for NGC 2768. The column descrip- 
tions are the same as in Fig. lAll Because a mask was observed 
twice independently on separate nights, some positions are as- 
signed two separate measured values. 

CaT value reached beyond the range covered by the V03 SSP models. 
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Table A3. Individual values for NGC 4494. The column descrip- 
tions are the same as in Fig. lAll 

*CaT value reached beyond the range covered by the V03 SSP models. 



